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ABSTRACT 



(N 

Context. The Milky Way bulge is the nearest galactic bulge and the best laboratory for studies of stellar populations in spheroids based 
on individual stellar abundances and kinematics. The observed properties point to a very complex nature, which is hard to extrapolate 
~"> ' from a few fields. 

rn \ Aims. We present a method to obtain reddening maps and to trace structure and metallicity gradients of the bulge using data from 

the recently started ESO public survey Vista Variables in the Via Lactea (VVV). The method is used to derive properties of the fields 

i-C ■ along the minor axis. 

ft' Methods. We derive the mean J - K s color of the red clump (RC) giants in 1835 subfields in the Bulge region with -8° < b < -0.4° 

and 0.2° < / < 1.7°, and compare it to the color of RC stars in Baade's Window for which we adopt E(B - V) = 0.55. This allows us 
to derive the reddening map on a small enough scale to minimize the problems arising from differential extinction. The dereddened 
magnitudes are then used to build the bulge luminosity function in regions of ~ 0.4° x 0.4° to obtain the mean RC magnitudes. These 
are used as distance indicator in order to trace the bulge structure. Finally, for each subfield the derived distance and extinction values 
have been used to obtain photometric metallicities through interpolation of red giant branch colors on a set of empirical ridge lines. 
The photometric metallicity distributions are compared to metallicity distributions obtained from high resolution spectroscopy in the 
same regions. 
^- ' Results. The reddening determination is sensitive to small scale variations which are clearly visible in our maps. Our results are in 

\& | agreement within the errors with literature values based on different methods, although our maps have much higher resolution and 

more complete coverage. The luminosity function clearly shows the double RC recently discovered in 2MASS and OGLE III datasets, 
hence allowing to trace the X-shape morphology of the bulge. Finally, the mean of the derived photometric metallicity distributions 

l/~) ' are in remarkable agreement with those obtained from spectroscopy. 

Conclusions. The VVV survey presents a unique tool to map the bulge properties by means of the consistent method presented here. 
The remarkable agreement between our results and those presented in literature for the minor axis allows us to safely extend our 
method to the whole region covered by the survey. 
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rN 1 ■ Introduction for a high complexity of the bulge structure and stellar popula- 

3 , tion cont ent: i) The Milky Way bulge has a radia l metallicity 

Several studies have attempted to address the problem of the g^gm dZoccali et all l2"008l I.Tohnson et all I20TI most likely 

origin of the Galactic bulge, trying to reconstruct its star for- flatten ing in the inner regions dRJche't al. 2007); ii) it seems 



mation history through the analysis of b oth stellar age s and 
chemical abundances (McWilliam & Rich 1994; Zoccali et al 
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Fulbright et al. 2007; Rich et al. 2007; Lecureuretal 
2007; iZoccali et alj 120081; Idarkson et alj 120081; iBrown et al 

"ID 



to host two kinematica l ly and chemica l ly dis tinct components 
dBabusiaux et all 120 10: Gonzale z~eTaT1 l201lh: and hi) it has 
an X-shape (McWi lliam"& Zoccalill2010l;lNataf et al.1l201f3) . All 
these observations strongly argue against the assumption that the 



20101 Uohnson et alj EQUJt [Gonzalez etalj |2QJJ. Early stud- p r0 p er ties of the Galactic bulge stellar population in a few low 
les were restricted to a low extinction window at a latitude extinction windows could be extrapolated to the whole bulge. 
b = -4° along the minor axis of the bulge, the so-called Baade's Further ana i ysis of the chemical properties of stars in other bulge 
Window. However several recent papers presented the evidence regionSj away from the minor axiSi has been partia i ly obtained 
from the latest sample of microlensed dwarfs from Bensb vet al.l 



Send offprint requests to: O. A. Gonzalez (120101) . Although microlensing provides chemical abundances 

* Based on observations taken within the ESO VISTA Public Survey all across the bulge, a larger sample is required in order to verify 
VVV, Programme ID 179.B-2002 
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the presence and extent of gradients. Consequently, at present, a 
different approach based on a larger coverage with a statistically 
significant stellar sample is the key to disentangle the general 
properties of the bulge. 

In this context, a base for such larger coverage studies is 
required i n order to ana l yze a nd correctl y interpret results. In 
particular, IZoccali et al.l (120031) and later iJohnson et all {2011) 
showed how the upper red giant branch photometry can be used 
to obtain photometric metallicity distributions for a large num- 
ber of bulge stars in two fields at b — -6° and b - -8°, re- 
spectively. Using a similar photometric approach to derive the 
metallicity distribution of a much larger bulge region will pro- 
vide for the first time a general picture of the iron content of 
the galactic bulge. Unfortunately, in order to reach this goal a 
first challenge arises from the extinction variations across dif- 
ferent bulge sightlines. While most of the studies re lied on the 
large scale extinction maps of Schleg el et alJ (11998), it is well 
known that these maps become unreliable when inner regions 
of th e bulge closer to the g alactic plane (\b\ < 3) are consid- 
ered. Schultheis et al. ( 1999), based on DENIS photometry, and 
Dutra et al. (2003), based on 2MASS, used the red giant branch 
photometric propert ies to obtain 2D extinction maps in several 
inner bulge regions. Marshall et al.1 (120061) presented 3D models 
of the extinction properties of the inner Galaxy with a resolution 
of 15'. However, at this point a homogeneous extinction map of 
the bulge including both the outer and the most inner regions, 
where high resolution is required to avoid differential extinction, 
is still necessary. 

Re garding the str ucture of the bulge, t he work of iNataf et alJ 
d2010 ) and Mc Will iam & Zoccalil d2010l) , based on OGLE and 
2MASS photometry respectively, provided hints for a X-shape 
morphology in the outer regions of the bulge which however, 
seems to merge into a single structure in the innermost regions, 
i.e. the bar (Rattenburv et al. 2007). Both Mc William & Zoccali 
(2010) and Nataf etal. (2010) used the observed properties of 
the RC stars at different lines of sight, which unfortunately re- 
stricted the studies to the outer bulge regions due to the non- 
complete optical coverage of OGLE and the relatively bright 
limit magnitude of 2MASS. For a thorough study of the bulge 
structure, and to homogeneously trace both the inner and outer 
regions, high spatial coverage and deep near-IR photometry are 
necessary. 

Both reddening and structural properties can be obtained 
from the photometric analysis of RC giants across the Galactic 
bulge. RC giants are core-helium burning stars, equivalent 
of the horizontal branch stars in old metal-poor populations, 
and are easily identified in bulge color-magnitude diagrams. 
Their intrinsic magnitude has a fairly well understood de- 
pendence on population properties such as age or metallicity 
dGirardi & Salarisl 1200 ll) . and reddening effects in color can 
be clearly identified. For old and relatively met al-rich popu- 
lation such as found in the Galactic bulge (e.g. iBrown et alJ 
2010, and references therein) the population corrections for red 
clump absolute K-band magnitude are of the order of 0. 1 mag 
(Sala ris & Girardil2002l) . For this reason, photometric properties 
of red clump giants can be used to trace simultaneously extinc- 
tion and relative distances in the bulge. The VISTA Variables in 
the Via Lactea (VVV) Survey, mapping 52 deg 2 including th e 
southern disk and bulge of the Milky Way (Min niti et al.ll2010l) . 
provides for the first time deep enough near-IR photometry to 
trace RC stars in all bulge regions including highly reddened in- 
nermost parts. 

In this paper, we present the first analysis of the Galactic 
bulge RC giants based on VVV data for a region along the mi- 



Table 1. Description of the analyzed fields. 



VVV 


/ 


b 


Obs. 


seeing (arcsec) 


field 


deg 


deg 


Date 


J 


H 


K, 


b236 


0.93 


-7.60 


13/08/2010 


0.90 


0.80 


0.79 


b250 


0.91 


-6.51 


13/08/2010 


0.93 


0.83 


0.88 


b264 


0.89 


-5.41 


14/08/2010 


0.91 


0.83 


0.81 


b278 


0.88 


-4.32 


22/04/2010 


0.91 


0.86 


0.84 


b292 


0.87 


-3.23 


30/08/2010 


1.20 


1.09 


1.13 


b306 


0.86 


-2.14 


30/08/2010 


1.41 


1.10 


1.14 


b320 


0.85 


-1.04 


22/04/2010 


0.81 


0.76 


0.71 



nor axis. We analyze the RC properties in the color magnitude 
diagrams along different lines of sight to infer reddening prop- 
erties and mean distances. This information is then coupled with 
2MASS photometry of the upper red giant branch (RGB) stars in 
order to further derive bulge metallicity distributions. This paper 
describes in detail the method and compares the results with the 
present literature. In a forthcoming paper we aim to apply the 
method to the complete 300 deg 2 coverage of the VVV bulge 
area, deriving homogeneous extinction, morphological structure 
and metallicity maps of the entire Galactic bulge. 

2. VVV survey data 

The VISTA Variables in Via Lactea (VVVfl is the Galactic 
near-IR ESO VISTA public survey that started collecting data 
in 2010. The VVV area covers ~ 520 deg 2 and will be imaged 
repeatedly in /f r band with the aim to search for variable stars 
such as RR Lyrae in the Milky Way bulge, as well as Cepheids 
and eclipsing binaries across the plane. In addition to the main 
variability campaign, the survey area is fully imaged in 5 photo- 
metric bands: Z, Y, J, H, and K s , 

The VVV survey has two main components: the 300 sq. deg 
VVV Galactic bulge survey area covering -10° < / < +10° and 
-10° < b < 5°, and the 220 sq. deg VVV plane survey area 
spanning between 295° < / < 350° and -2° < b < 2°. In the first 
observing season in 2010 95% of the total area has been imaged 
in J, H and K, wavebands, while the YZ imaging covered 84% 
of the disk and 40% of the bulge areas. The variability campaign 
has also started, but mostly it will be carried out in the upcoming 
observing seasons. The survey is expected to span 5 years of 
observations, enabling to characterize well both short and long 
period variable stars in the Milky Way bulge and the plane. 

Here we describe only data used in the analysis presented 
in this paper: the single-epoch JHKs coverage of the bulge. The 
more complete description o f the whole survey a nd the observa- 
tion strategy is available in Minn iti et alJ (1201 0[) and some firs t 
results from the VVV survey are presented in lSaito et alJ (12010). 

The observations are carried out in service mode with the 
VIRCAM camera on VISTA 4.1m telescope (Emerson et al.l 
12004 lEmerson & Sutherland! 1201(H) located at ESO Cerro 
Paranal Observatory in Chile on its own peak, some 1 .5km away 
from the Very Large Telescope. The bulge survey area consists of 
196 tiles, each 1.48 x 1.18 sq. deg in size. VIRCAM focal plane 
is sparsely populated with 4x4 Raytheon VIRGO 2048 x 2048 
pixel arrays, and 6 exposures (so called pawprints) are needed to 
fully cover the 1.5 square degree field (so called tile). The mean 
pixel scale is 0.34 arcsec, and the point spread function deliv- 
ered by the telescope and camera has FWHM ~ 0.51 arcsec. 



http://vvvsurvey.org 
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Fig. 1. Distributions of the RA and DEC differences in arcsec 
between matched sources of single band, J and Ks, photomet- 
ric catalogs. Crossmatching for field b278 is shown here as an 
example 



2 |- 

2 §- 

4 '- , 
6=- ( J 



-J„J = 0.029 




13 
Ks„. 



Fig. 2. J,H and K s magnitudes differences between VVV cat- 
alogs and 2MASS. Dashed lines show the mean difference be- 
tween both catalogs in the range denoted by the solid lines 



The median image quality measured on the reduced VVV tile 
images ranges between 0.8" for K s , 0.9" for J band, and up to 
1 .0" for Z-band. Few tiles were observed out of seeing constrains 
and were later repeated. This is the case for tiles b292 and b306 
for which we currently only have reduced images taken with a 
seeing larger than 1 .0". In the future analysis of the whole bulge 
sample, these tiles will be replaced with newer observations. 

VVV tiles are made combining 12 (two per pawprint) images 
for H, K s and J band tiles. Each H and K s image is a single 4 sec 
exposure, while J band images are sums of two 6 sec exposures. 
Each pixel in a tile has therefore a minimum of 16 sec exposure 
in H and K s and 24 sec in /, and is an average of at least 4 indi- 
vidual exposures. VISTA data are processed through the VISTA 
Data Flow System (VDFS) by the Cambridge Astronomical 
Survey Unit (CASU). The basic data reduction steps consist of 
reset correction, dark subtraction, linearity and flat field correc- 
tion, sky background subtraction, destriping step, jitter stacking 
and tiling. Finally the astrometric and photometric calibration 
are applied using 2MASS sources in the images, and the single 
band catalogues are produced (Lewis et al. 2010). In this work 
we use the single band tile images and catalogues version 1 . 10. 

The tiles analysed in this paper are located along the minor 
axis at I ~ 0° and spanning -8° ^ b < -0.5°. Table[T]lists center 
tile galactic coordinates, date of observation, and the FWHM of 
the PSF in arcsec measured on tiles. Given that data for each tile 
in three filters J, H, and K s are taken consecutively, the airmass 
is approximately the same for all filters in a tile. 

2. 1 . Multiband catalogs 

Target coordinates and J,H,K S magnitudes were extracted from 
the single band catalogs produced by the CASU. These cata- 
logs were then matched using the code STILTS (Tavlor 2006). 
In particular, the routine tskymatch2 is used to perform a match 
based on the proximity of sky positions of stars in each cata- 
log. Several tests were carried out by modifying the maximum 
separation radius to match detections between J, H and K s cat- 



alogs. The majority of stars were always found within a separa- 
tion smaller than 0.5", in good agreement with the astrometric 
precision expected from CASU pipeline. Therefore, a maximum 
allowed separation of 1" was used to ensure that our matched 
catalogs contain most of the sources without risk of including 
wrong matches at a larger separation, and at the same time we 
still leave the possibility for some minor astrometry distortions 
in the large VIRCAM field of view. Figure Q] shows the distri- 
butions of the RA and DEC difference for J and Ks catalogue 
matches in the VVV field b278 . 

The single band tile catalogs produced by CASU contain a 
flag which identifies detections as stellar or non-stellar sources. 
We decided to only use sources with the stellar flag in all three 
bands. Although this cleaning procedure implies that only ~75% 
of the original number of detections remain in the final multi- 
band catalog, it allows us to work with as clean as possible sam- 
ple of bulge stars, which is ideal for our purpose, and it does not 
bias our results. 

2.2. Photometric calibration 

The VVV images extend about 4 mag deeper than 2MASS0 in 
spite of the short exposure times because of the difference in 
telescope apertures and spatial resolution. However, the larger 
telescope aperture of VISTA results in saturation for brightest 
stars. For this reason 2MASS photometry is required to com- 
plete the bright end of the VVV catalogs. Figure [2] shows the 
J, H and K s magnitude differences between VVV and 2MASS 
sources in a common 10x10 arcmin window within the tile b278. 
We only adopted 2MASS sources with best photometric quality 
(flag AAA). Two important points can be observed from Fig. [2] 
i) The saturation magnitude in the VVV photometry is clearly 
observed at K SYVV < 12; and ii) there is a very small offset be- 
tween the zero points of both catalogs. Still, we must ensure that 
both catalogues are in the same photometric system to be able to 
complete the saturated bright end of the VVV photometry with 
2MASS. In doing this, only stars common to both surveys and 
with 13 > /Ts.vvv > 12 (see Fig. 2) were used to bring the VVV 



2 More details about VISTA data processing is available from 3 The exact difference between VVV and 2MASS photometry depth 
http://casu.ast.cam.ac.uk/surveys-projects/vista/data-processinglepends strongly on the stellar density in the field. 
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Fig. 3. VVV observed (J - K S ,K S ) CMD of a 20' x 20' region 
within the tile b278 calibrated and completed with 2MASS pho- 
tometry. 



catalogues onto the 2MASS photometric system. The magnitude 
range selected for the calibration ensures that the differences are 
not affected by larger errors in the photometry for 2MASS or sat- 
uration for VVV measurements. This procedure was carried on 
independently for each VVV tile and zero points were calculated 
for each one of the catalogs. Once calibrated, stars brighter than 
K s ~ 12 from 2MASS were used to complete the bright end of 
our final VVV catalog. These catalogs, corrected for saturation 
and calibrated to 2MASS zero points were used in further anal- 
ysis. Hence, hereafter we refer to the usual J,H,K S magnitudes 
corresponding to the 2MASS photometric system. 

2.3. Bulge minor axis color-magnitude diagram 

As an example, Figure [3] shows the K s , J - K s color-magnitude 
diagram (CMD) of a selected region in the field b278. The pho- 
tometry samples the entire RGB from the base up to the tip 
(K s ~ 8). The RGB is quite bent as expected for metal-rich 
populatio ns and its relative ly large width is due to metallicity 
spread dZoccali et al.ll2.003l) . depth effects and some differential 
reddening. The red clump is located at K s ~ 13, and it is par- 
tially merged with the RGB bump. Although the photometry is 
deep enough to properly sample the sub giant branch, the almost 
vertical bluer sequence due to foreground disk stars prevents to 
detect its exact location. 



3. Relative reddening values 

The extinction corrections are derived by using a differential 
method based on the comparison between the mean color of the 
RC stars observed in all fields with those of a reference field with 
well known extinction, the Baade's Window field. The important 
assumption we are taking here is that the differences in terms of 
average age and metallicity of the bulge stellar populations can 
be accounted for with only a small correction on the color of the 
RC stars in the bulge fields, given the old age and nearly solar 
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Fig. 4. Right panel shows the observed color magnitude diagram 
for a 10 x 10 arcmin region centered in Baade's Window and a 
selection box corresponding to the red clump. Left panel shows 
the color distribution of these red clump stars. Overplotted is the 
best fit Gaussian curve with the mean RC color. 
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Fig. 5. Red clump region of the observed color magnitude di- 
agram for tile b306 (right panel) compared to the de-reddened 
one obtained by our procedure (left panel). Only a subsample of 
the stars are plotted in each CMD and contour levels are shown 
to mark the clump position and shape. The arrow in the right 
panel shows the reddening vector for a change in E(B-V) of 0.5 
magnitudes. 



metallicit y of the Bulge and the presence of a known metallicity 



metallicity or tne tsulge ana tr 
gradient dZoccali et al.l l2008). 



In order to obtain the properties of a reference RC popula- 
tion we have analyzed one of the fields for which we have spec- 
troscopically determined iron and alpha element abundances 
(Zoc cali et alj|2008c iGonzalez et alj|201 ll) . We have measured 
the color of the red clump in a 10' x 10' region centered in 
/ = 1.14 and b - -4.18 (hereafter Baade's Window) for 
w hich we have adopte d a reddening value of E(B-V)=0.55 as 
in IZoccali et al.l (120081) . Figure [4] shows the observed CMD for 
this region and the (J - K s ) color histogram for which we find 
a mean value of (/ - K s )™=0.96 based on a Gaussian fit to the 
color distribution. Adopting this same procedure, the extinction 



Gonzalez et al.: Reddening and metallicity maps of the Milky Way Bulge 



E(B-V) of a bulge field at (l,b) can be related to that of Baade's 
Window with the following equation: 



E(B - V) = E(B - V) BW - A(J - K s ) RC /(0.fn - 0.35) 



(1) 



where E(B- V)bw is the extinction in Baade's window and A(7- 
K s )rc is the difference between the color of the RC in Baade's 
Window, (J - K s )^, and the observed color of the RC for the 
field at (Lb), (J - K,)'k r , as obtained from the Gaussian fit. For 
our pr ocedure we adopted the reddening law from lCardelli et alj 
(fl989h : A K = 0.35 ■ E(B - V), A H = 0.59 • E(B - V) and A, = 
0.87 ■ E(B - V). We note that different reddenin g laws have been 
deriv ed for some inner bulge fields (\b\ < 2) (Nishivama et al.l 
2009, and references therein), however these do not cover the 
full area that we are investigating. On the other hand, the bulge 
RR Lyrae study from lKunder et al.l (120081) showed that, although 
variations are observed in particular directions, a standard red- 
dening law of Ry ~ 3.1 is, on average, valid for bulge studies 
at larger distances from the Galactic plane. Sti ll, the adoption 
of a di fferent reddening law, such as the one in Nishivama et al.l 
(2009) where A K = 6.528 ■ E(J - K) instead of the standard 
A K - 0.640 ■ E(J - K) would affect our results by less than 0.06 
mags in K s band, for the range of extinctions observed in the 
regions where we provide a bulge stellar population analysis. 

We are able to use this technique in all bulge fields cov- 
ered by the survey as VVV photometry reaches the RC even in 
the highly extincted regions. We have divided each VVV field 
into smaller subfields of 10' x 10' for (b < -5°), 6' x 6' for 
(-5° < b < -1.5°) and 3' x 3' for (-1.5° < b < 0°) in order 
to minimize the effects of differential reddening and at the same 
time, to accurately trace the RC with robust number statistics. 
For fields at larger distances from the plane (b < -5°) the den- 
sity of stars and the differential reddening variations are smaller 
while the opposite occurs in regions closer to the plane. For this 
reason we use sub-fields with slightly different sizes and derive 
a single reddening value for each of the sub-fields. 

A fully de-reddened CMD zoomed on the RC region is 
shown in Fig [5] for tile b306. The effect of differential redden- 
ing is clearly observable as an increase of the RC width (Fig. [5] 
right panel), and it is consistent with the reddening vector shown 
in the figure to indicate the effect of E(B-V) difference of 0.5 
mag. We have carried out the same procedure in all 6 tiles along 
the bulge minor axis (Table[T]l. The main objective was to obtain 
a homogeneously de-reddened dataset of J,H,K S photometry to 
investigate RC properties and metallicities, however at the same 
time we also got a consistent reddening map for this bulge re- 
gion. 



3.1. Minor axis reddening map 

Given the density of stars in the inner fields along the minor axis, 
the area was split into 1835 subfields with boxes of 3' on a side 
for fields within inner 1.5°, 6' for -5 < b < -1.5° and 10' for 
b < -5°. Our map is therefore in-sensitive to reddening varia- 
tions on smaller scales. Figure [6] shows the reddening map ob- 
tained for the bulge minor axis between b - -8° and b - -1.5° 
at longitudes between 0° and 2°. As expected, although the small 
scale variations are clearly observed, they become particularly 
strong in the inner regions as can be observed in Figure|7]which 
shows the reddening map for b > - 1 .4. This indicates that the 
use of lower resolution maps such as the one of Marshall et al. 
(2006) (15') for the very inner bulge regions, near to the galactic 
plane, will not be sensitive to such variations. 
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Fig. 6. Reddening map along the minor axis of the bulge. 
Regions with b < -5° are smoothed to 6' boxes for visualiza- 
tion. 
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Fig. 7. Reddening map for the inner regions of the bulge along 
the minor axis. 



High resolution extinction maps for the g alactic bulge, in- 
cludin g the i nner regions , have b e en obtained by Schultheis et al.l 
(199 9I) andlDutraetalJ (l2003h . ISchulfheis et al.l (Il999l) used 
Bert elli et alj (U994) isochrones fitting to red giant branch stars 
in DENIS (J,H,K S ) CMDs built in subfields of 2 x 2 arcmin. 
Their reference isochrones where calculated by assuming a 10 



Gyr solar metall icity bulge populat ion at 8 kpc from the Sun. On 
the other hand, iDutra et all (120031) determined reddening with 
a similar procedure based on 2MASS photometry, but using a 
composite observed bulge CMD as reference instead of adopt- 
ing an isochrone model. 
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Fig. 8. Compa rison between th e A<- va lues obtained in our work 
and those o f Schulthei s et al.l d 19991 black filled circles) and 
iDutra et all d2003l red crosses) 



Figure [8] sh ows the differenc e betw een o ur reddening value s 
and those from lSchulfheis et alj d 19991) and IDutra etaIT(l2003l) . 
The comparison was made by matching the central coordi- 
nates of the s ubfields of the different studies. We have applied 
ICardelli et alj (1989) extinction transformations to calculate Ak 
values from these studies, in which other transformations were 
adopted, in order to obtain a consistent comparison with our 
results. We systematically find larger reddening values and the 
difference becomes particularly conspicuous in high extinction 
regions. This can be caused by the different assumptions that 
each technique adopted. Our map is relative to the zero point 
adopted from the observed (J-K s ) color of the RC in a 10' x 10' 
field at Baade's Window, for which we assume a reddening 
value of E(B-V)=0.55. This alone could introduce a system- 
atic offset depending on this particular reddening value. Also, 
our method assumes that the intrinsic color of the RC is, in the 
mean, the same all across the Bulge. Unfortunately, this or sim- 
ilar assumptions are always n ecessary. In particu lar, the studies 
of lSchultheis et all d 19991) and lDutra et~al] d2003l) are also based 
on assumptions of the homogeneity of stellar populations in the 
bulge and refer to particular mean values in age and metallicity 
which might introduce the observed offsets. However, RC prop- 
erties used in our work have a smaller dependence on metal- 
licity, within the known bulge values, in comparisson to those 
of the RGB. Altogether, when considering the difference in the 
methods and datasets, the agreement between the three studies 
is remarkable up to Ak ~ 1 . For higher extinctions, given that 
those regions are also the ones suffering higher crowding prob- 
lems, we expect that the high resolution and deeper photometry 
of VISTA provide the most reliable results. 



4. Bulge luminosity function and distances 

With the de-reddened photometry of VVV obtained as explained 
in the previous section, we are abl e to use the RC magnitu de as a 
distance indicator (P aczynski & Staneki ll998l:lAlves et a"i]|2002t 
iGrocholski &S araiedini 20021) We divided each VVV field in 
four subfields of ~ 0.4° x 0.4° in which we built the luminosity 
function (LF) and used it to obtain the mean observed magnitude 
of the RC, ^f c , by means of a Gaussian fitting procedure (e.g. 




Fig. 9. Lower panels show the luminosity functions for a subfield 
in Baade's Window and at b=-6. The underlying RGB is fitted 
with a second order polynomial and the red clump is fitted with 
two Gaussians centered at K„ - 12.9 and K Sg = 13.2 for b=-6 
and with a single Gaussian centered at K Sa = 12.9 for Baade's 
Window. Upper panels show the corresponding CMD for each 
field, oriented with color along the y-axis and magnitude along 
the x-axis. Density contours denote the single and double RC in 
each field. 



IStanek & Garnavichlll998l;lBabusiaux & Gilmorell2005l) follow- 
ing the function: 



N(K S0 ) = a + bK S0 



cKl + 

^0 



N, 



RC 



O-RC 



V2^ 



exp 



r« c 



^o) 2 



2ai 



(2) 



We first applied a color cut at (J - K s ) > 0.5 to minimize the 
contamination from foreground stars and then fitted the under- 
lying RGB LF using a second order polinomial corresponding 
to the first 3 terms of Eq. (2). The RC was then fitted with two 
Gaussians, in the form of the forth term in Eq. (2), in the regions 
where the two red clumps are observed, i.e for latitudes b < -5° 
(Mc Willia m & Zoccalil l2010) and with a single Gaussian in the 
inner regions. An additional Gaussian fit was performed when 
the red giant branch bump (RGBb) was detected. Upper panels 
of Fig|9]show the final de-reddened CMDs for Baade's Window 
and the field at b — -6° where the density contours show clearly 
the single and double RC features in each field. The lower pan- 
els of Fig. [9] show the de-reddened K So band luminosity func- 
tions for both fields with overplotted polynomial + Gaussian 
fits. The polynomial function fits very well the underlying red 
giant branch population, while Gaussian functions are required 
to fit RCs and RGBb. A single RC is observed at K S(t -12.9 
in b - -4°, while at b - -6° the RC region shows two com- 



ponents at K So ~ 13.2 and K S( 



is observed at K s , 



13.6 in/? 



12.6. Additionally, the RGBb 
4° . The lat ter is in agreemen t 



with the measurements of the bulge RGBb bv lNataf et al.l d201ll) 
using OGLE photometry. 

We use the observed mean RC magnitudes in order to ob- 
tain distances for each of our subfields through comparison of 
the measured values with the intrinsic RC magnitude for a 10 
Gyr solar metall i city p opulation of Mk s = -1.55 based on 
Pietrinferni et al. (2004) Teramo stellar evolution models (see 
also Ivan Helshoecht & Groenewegen 2007). 

As we already discussed in the case of the reddening cal- 
culations, the assumption of homogeneous stellar populations 
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on the bulge are the main source of uncertainties. The intrin- 
sic magnitude of the RC is known to have a dependence on age 
and metallicity (Salaris & Girardi 2002; Pietrzyhski et al. 2003). 
While the age of the bulge is expected to be homogeneously 
old (~10 Gyr), the metallicity is kn own to change at least along 
the mi nor axis (Zoccali et al. 2008). According to IZoccali et al.l 

(2008), the Bulge metal content varies from [Fe/H] 0.3 at 

b = -12° to [Fe/H] 0.1 at b = -4° a nd there is a hin t 

for a lack of gradients on the inner regions dRichetal.1 12007). 
Teramo stellar evolution models indicate that a 0.2 dex variation 
in metallicity would affect the intrinsic mag nitude of t he R C by 
~0.10 mag, consistent as well with Salaris & Girardi (2002). If 
for a solar metallicity population at 8 kpc from the Sun we as- 
sume a main metallicity 0.20 dex lower, our method will result 
in measurement of distance of ~8.2 kpc instead. Therefore, al- 
though this introduces a small error into our calculations we are 
still able to trace the large scale structure of the bulge. 

iMc William & Zoccalil d2010l) and lNataf et all d2010l) showed 
for the first ti me the presence of a double RC observed at 
large latitudes. Mc William & Zoccali d2010f) argued that this 
feature could be consequence of a X-shaped Bulge. Furthermore, 
Saito et all d201ll) mapped the density of RC stars across the 
bulge using 2MASS photometry and concluded that the X-shape 
morphology merges to the bar in the inner regions. Our results 
are in agreement with the measurements in previous studies, 
showing the clear double peak in the LF of the bulge as can 
be seen in Fig(9] Our method can be then used to map the in- 
nermost Bulge re gions which are not considered in the work of 
Saito et all d2011l) because of the limiting depth and blending in 
2MASS photometry in these inner regions. This analysis is un- 
derway and will be the subject of a dedicated paper. 

5. Bulge photometric metallicity distributions 

The final goal of our work is to use the information obtained in 
the previous sections such as reddening and mean distance to 
study stellar population metallicity gradients in the bulge. We 
compare the derived field CMDs with an empirical grid of RGB 
ridge lines of globular clusters in the [M^ s , (J - K s )o] absolute 
plane. The empirical templates of Galactic halo and bulge glob- 
ular clusters (M92, M55 , NGC 6752, NGC 362, M 69, NGC 
6440, NGC 6528 and NGC 6 791), with well kno w metallicity, is 
selected from the sample of Valenti et al. (2004). We then derive 
a metallicity value for each field RGB star from its (J—K s )q color 
interpolating within the empirical grid. This proce dure has been 
prove n to be effe ctive in b ulge globular clusters (Valenti et al. 
120071) . field dZoccali et all [2003: Joh nson et all 120 111) as well 
as in studies of mo re distant galaxies (e.g. Rejkuba et al. 2005; 
lLianouetal]l201ll) . 

To transform the observed field CMDs into the absolute 
plane we use the distance and reddening values derived in the 
previous sections, while for t he cluster ridge fine s we adopt the 
estimates listed in Table 3 of IValenti et al.l d2004h . It is worth to 
mention that the cluster templates are derived from a homoge- 
neous photometric datab ase calibrate d onto the 2MASS photo- 
metric system (see lValenti et alJl2004 for more details) and they 
have been selected in order to span a large metallicity range, be- 
tween -2.16 and +0.35. In particular, we have included in our 
empirical grid the RGB ridg e line of NGC679 1 which is one 
of the oldest (age 6-12Gvr. ICarnev et alj 2005|, and refere nces 
therein) and metal-rich ([Fe/H] =0.35, (Oridia et al.ll2006l and 
references therein) open cluster. Furthermore, we notice that the 
alpha-abundance patterns of the selected templates - enhanced 
for metallicity up to that of NGC6528 and solar for higher metal- 



licities - match very well the abundance patterns observed in 
the Bulge fields (Gonzalez et al.ll201 lb lAlves-Brito et~ai1l2010l: 
Lecureuret al. 2007; Fulbright et al. 2007). This implies that 
there is no need to adjust the fiducial ridge line of the clusters 
to match the alpha-abundances of the fields. In the interpolation, 
stars brighter than M^ = -4.5 were used to include only the 
RGB region most sensitive to metallicity variations, and to avoid 
the use of red clump stars (See Zoccali et al. 2003, for a detailed 
discussion). Additionally, a color cut was applied including only 
stars with (J - K s )q > 0.6 to minimize the disk contamination. 

Although this method has been shown to work well in spe- 
cific cases (i.e. single stellar population and distant galaxies), 
here a problem may arise from the assumed distance modulus 
along different latitudes. First, as discussed in the previous sec- 
tions, the distance estimates are derived assuming an intrinsic 
RC magnitude which depends on the metallicity. However, the 
expected error of 0.10 mag for RC stars due to metallicity gra- 
dient has small effect on our photometric metallicity measure- 
ments because these are based on RGB morphology which has 
significantly larger sensitivity to metallicity variations. To show 
that we compare in the next section our photometrically deter- 
mined metallicity distributions with iron abundances from high 
resolution spectroscopy. 

A second and more important problem arises from the ob- 
served double red clump in the luminosity function, which traces 
the X shape of the bulge. The choice of the bright or the faint RC 
as the observed value to obtain the absolute magnitudes has a 
strong impact on the results because it implies a change of up to 
1 .0 mag along the minor axis. Clearly, whichever RC we assume 
to be the dominating population, RGB stars from the popula- 
tion corresponding the other RC will have an error in magnitude 
as large as the separation of both clumps. Since it is impossi- 
ble to distinguish the two populations on the RGB, we adopt 
the weighted average magnitude of the two RCs as the observed 
magnitude to calculate the distance modulus for the RGB stars 
in that field. This magnitude was calculated as: 



Kf = (Kf o cl N RC i + Kf o c2 N RC 2)/(N RCl + N RC2 ) 



(3) 



where N R ci and N R ci are the number of stars in each clump with 
mean magnitudes ^f cl and K^ 2 respectively. 

Therefore, the errors for metallicities for individual stars 
obtained using our interpolation method are expected to be 
significantly larger than those based on high resolution spec- 
troscopy. However, the photometric metallicity method provides 
two unique benefits: (i) the large number statistics, and (ii) most 
importantly the large coverage, which allows us to investigate 
gradients. 

5.1. Comparison to spectroscopic MDF and the minor axis 
gradient 

The implications of the errors in the method due to assump- 
tions on the reddening law, on the selection of stars on the RC 
and the assumptions on population corrections for RC absolute 
magnitude, are hard to take into account within our calculations 
rigorously and therefore a direct comparison with a more pre- 
cise method such as high resolution spectroscopy is particularly 
useful. Several spectroscopic studies have been carried out in 
fields along the minor a xis of the bu l ge. In particular the metal- 
licity distributions of Z occali et alJ ^2008) presented the evi- 
dence for metallicit y gradient between la titudes -4°, -6° and 
-12°. The study of iJohnson et al.1 (1201 ll) provided the metal- 
licity distribution of a field at b = -8°, which is consistent 
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Fig. 10. Left panels: Comparison of the photometric (hashed histogram) and spectroscopic (red dotted histogram) MDFs along the 
Bulge minor axis at b — -4°, -6° and -8°. The average [Fe/H] and dispersion cr s are also shown for all MDFs. Right panels: 
CMDs in the absolute plane of the fields located along the minor axis at b — -4°, -6° and -8° with overplotted the empirical RGB 
templates. Thick black circles refer to the stars used to derive the photometric MDFs. 



wi th the expectat i on fro m the metallicities and gradient observed 
by IZoccali et alj d2008l) . To check the reliability of our photo- 
metric metallicity distribution functions (MDFs), we compare 
the MDFs from our photometric metho d in ~ 30' X 30' fields 
with those from spectroscopi c studies of | Zoccali et all d2008l) at 
b = -4° and b = -6° and of I Johnson et all d201 lh atT= -8°. 

Figure[l0]shows the remarkable agreement for the mean pho- 
tometric (hashed histograms) and spectroscopic (red dotted his- 
tograms) MDFs. The metallicity gradient is clearly observed in 
our data, and the shape of the photometric distributions resem- 
ble very well those from spectroscopic studies. The dispersion 
values in the photometric [Fe/H] distributions are larger than 
the spectroscopic ones. This may be produced by a combina- 
tion of factors. In particular, our method relies on extrapola- 
tion for stars with metallicities larger than 0.35, the value of 
the most metal-rich template corresponding to the open cluster 
NGC6791. This could generate a metal rich tail with no reli- 
able abundances and which might increase the dispersion. On the 
other hand, the metal-poor tail has a possible disk contamination 
which could end up also increasing the dispersion of our distri- 
butions. However, we do not exclude that such metal-poor stars 
could actually belong to the bulge. Altogether, we demonstrate 
that our metallicity determinations are as reliable in the mean as 
the high resolution spectroscopic studies. This presents a power- 



ful tool to trace how metallicity varies not only in discrete fields 
along the minor axis, but all across the outer bulge. Figure QT| 
shows a map of mean metallicity in bins of ~ 0.4° x 0.4° and 
it confirms the metallicty gradient along the minor axis. These 
maps are restricted to \b\ > 3° as the upper RGB in our dataset 
corresponds to 2MASS photometry which is limited by crowd- 
ing at inner regions. In the future paper we will present the metal- 
licity maps of the entire VVV bulge area for \b\ > 3°. 

6. Conclusions 

We used the VVV survey data to investigate the properties of 
the Galactic bulge along the minor axis. Within this bulge re- 
gion there have been several spectroscopic studies which we 
use to validate our method. The method is based on the use 
of properties of the red clump giants. The RC color is used to 
derive the reddening map along the minor axis extending from 
-8° < b < -0.5°. These regions suffer from high extinction and 
variations on very small scales are observed. Therefore a consis- 
tent bulge reddening map obtained from a homogenous tracer all 
across the bulge, from the innermost crowded and extinct regions 
to the outer parts is very valuable. Our results for th e extinction 
also sh ow good agreement with m aps published by iDutra et all 
(120031) . and lSchultheis et ai1dl999h . 
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Fig. 11. A map of the mean values of the metallicity distributions 
along the bulge minor axis for subfields of ~ 0.4° x 0.4°. 



We have used the de-reddened magnitudes to build the Ks 
band luminosity functions and to trace the RC along the minor 
axis. When using the RC as a distance indicator, we find excel- 
lent agreement with previous findings of a double red clump for 
latitudes larger than b < -5° which trace t he X shape of th e 
Galactic bulge dMc William & Zoccalill20Tol: iNataf et al.ll2010h . 
We also detect the additional bump in t he lum inosity function 
interpreted as the RGBb by INataf et all d2011l) . This study will 
be extended to the complete 300 deg 2 bulge region and will pro- 
vide a consistent trace of the structure of the bulge including the 
inner regions. 

Finally, we used the RC distances and reddening map to 
build the color magnitude diagrams in the absolute plane and 
then to obtain individual photometric metallicities for RGB stars 
in the bulge throug h interpola t ion o n the RGB ridge lines of clus- 
ter sample from Valenti et al. (2004). The final metallicity distri- 
butions w ere compared to tho se o btained in the spectroscopic 
studies of Zoc cali et al.l d2008l) and lJohnson et alj d2011l) show- 
ing a remarkable agreement. We demonstrated clearly the possi- 
bility to track the observed metallicity gradient along the minor 
axis. This provides, within the errors of the method, the opportu- 
nity to obtain clues for metallicity gradients along other regions 
of the bulge. Reddening and metallicity maps for the VVV cov- 
erage will be released electronically in the survey website. 
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